Journal of Structural Biology

Cryo-electron tomography (cryo-ET) enables near-native visualization of subcellular architectures, yet applying it to moderately thick, multilayered tissues such as the retina is hampered by inadequate vitrification and inaccurate depth-targeting. Here, we developed PLCT, an integrated approach combining modified high-pressure freezing, cryo-ultramicrotome trimming, and plasma-based cryo-FIB milling to overcome these barriers. PLCT reliably vitrified <100 m retinal strips with minimal ice artifacts, navigates precisely to the outer plexiform layer using morphological landmarks, and produces high-quality lamellae suitable for high-resolution cryo-ET. Subtomogram averaging (STA) analysis identified microtubules at 16.33 [A] within retinal horizontal cell processes. Importantly, STA also resolved a 10-nm-diameter filamentous structure at 24.81 [A] in the same processes, featuring six peripheral strands surrounding an elongated central density with continuous intervening cavities, an architecture consistent with intermediate filaments. Together with its native localization and immunoreactivity, these features collectively identify the filaments as neurofilaments. Separately, 3D reconstruction of synaptic ribbons uncovered a previously unrecognized "mahjong tile"-like fine ultrastructure. These results demonstrate that PLCT-produced lamellae are of sufficient quality to support structural analysis in native tissue. Although demonstrated on retinal photoreceptor synapses as a proof-of-principle, PLCT is inherently generalizable, with its depth-navigation and vitrification strategies directly applicable to any multilayered tissues. This work establishes PLCT as a robust, reproducible platform for depth-resolved in situ cryo-ET of multilayered tissues.

This paper presents and critically reviews the results of original and some literature based experimental studies conducted by the authors last years on the structural organization of DNA in dormant (starvation stress), anabiotic dormant (4 HR treatment) E. coli cells, as well as the K12 {Delta}dps strain, which lacks the Dps protein (Dps null E. coli). The experimental data includes small-angle synchrotron radiation diffraction (SAXS) and transmission electron microscopy (TEM) data. Synchrotron radiation diffraction experiments on K12{Delta}dps cells allowed us to conclude that peaks at 44.3, 22.1, and 14.8 angstrom resolutions are associated exclusively with ordered DNA organization. Peaks at 44.3, 22.1, and 14.8 angstrom resolutions are also observed for samples of dormant (starvation stress) cells and anabiotically dormant cells. Therefore, this ordered DNA organization also applies to samples of dormant and anabiotically dormant cells. A model is proposed that considers the ordered DNA organization in the cell as a cholesteric liquid crystal. The powder diffraction pattern calculated based on this model is compared with experimental small angle X ray scattering (SAXS) data obtained on Dps-null cell samples. The model completely reproduces the key features of the experimental diffraction pattern from Dps-null cell samples. Accordingly, the cholesteric liquid crystal model corresponds to DNA packaging in dormant and anabiotically dormant cells. Cholesteric liquid crystal ordering should be further considered in all models of cellular DNA packaging. To address the question of which structural organization of DNA predominates in the cell: the cholesteric liquid crystal or nanocrystalline or whether they coexist and fully manifest themselves under different external conditions, it is necessary to utilize the latest methodological advances in structural analysis.

Ammonia oxidation initiates nitrification and is closely linked to microbial N2O production. Ammonia monooxygenase (AMO) catalyzes the first and rate-limiting step of nitrification and is widespread across evolutionarily distinct ammonia-oxidizing archaea (AOA) and bacteria (AOB). The ocean is the largest biome for AOA and AOB, which have distinct ecological niches and markedly different sensitivities to nitrification inhibitors. However, the lack of archaeal AMO structures and inhibitor-bound AMO complexes has hindered mechanistic understanding of the architectural, catalytic, and inhibitory divergence between these two enzyme systems. Here, we report high-resolution cryo-electron microscopy (cryo-EM) structures of marine archaeal AMO captured in active and inactivated states within its native membrane environment, together with inhibitor-bound structures of estuarine bacterial AMO. Archaeal AMO forms an unexpected cup-shaped homotrimer composed of eight subunits per protomer and exhibits substantial architectural divergence from bacterial AMO. Integrated structural, biochemical, kinetic, and computational analyses reveal distinct periplasmic architectures, copper-center organization, and hydrophobic channels between archaeal and bacterial AMOs for ammonium acquisition, catalysis and inhibitor response. These findings provide a structural and mechanistic framework for understanding how archaeal and bacterial AMOs have diverged to distinct ammonia-oxidizing strategies and inhibitor susceptibilities across environmentally important ammonia oxidizers.

Human cytochrome P450 2C9 (CYP2C9) is a hepatic microsomal enzyme involved in the oxidative metabolism of clinically important drugs, but the structural organization of its oligomeric assemblies outside crystallographic packing environments remains poorly understood. Here, we report the cryo-EM structure of human CYP2C9 determined under aqueous, membrane-free conditions at 3.31 Angstrom resolution. The structure reveals a C2-symmetric hexameric assembly organized as a dimer of trimers. Individual protomers retain the conserved P450 fold and heme-binding architecture observed in previously reported crystal structures, indicating that assembly formation does not substantially perturb the catalytic core. The hexamer is stabilized by defined intra-trimer interfaces involving the N-terminal region and residues around Trp212 and Phe482, together with inter-trimer interfaces involving Leu71 and the 220-227 loop. These interfaces are distinct from the crystal packing contacts observed in CYP2C9 crystal structures, demonstrating that the assembly is not a simple recapitulation of crystallographic packing. Notably, the inter-trimer interface is located near the FG-loop-containing surface previously implicated in membrane association. This suggests that the observed hexamer may represent a membrane-free association of two trimers through membrane-related surfaces, whereas the trimeric arrangement itself may be compatible with membrane-associated organization. The structure therefore provides a framework for investigating how trimer formation, membrane interaction and local conformational changes in the FG-loop region may influence CYP2C9 function.

Infectious bronchitis virus (IBV) is a member of the Gammacoronavirus genus responsible for respiratory illness and weakened eggshells in infected chickens, adversely impacting the poultry industry. Escaping innate immune detection during infection is crucial for coronavirus proliferation in the host. The production of double-stranded RNA during coronavirus replication triggers innate immune sensors to create an antiviral state within infected cells. To counter this response, coronaviruses employ nonstructural protein 15 (nsp15) endoribonuclease to degrade double-stranded RNA. Here, we use cryo-electron microscopy and biochemistry to characterize IBV nsp15 interactions with RNA. While the overall structure and active site of IBV nsp15 strongly resemble previous studies of nsp15 from other coronaviral genera, we note that double-stranded RNA contacts several non-conserved residues peripheral to the enzyme active site. Our data show that these residue positions can have strong impacts on RNA cleavage suggesting unique solutions for RNA engagement across coronavirus species. We also demonstrate a preference for IBV nsp15 to cleave double-stranded RNA over single-stranded RNA and observe nsp15 hexamers with two double-stranded RNAs bound simultaneously. This study reinforces the need to study diverse coronavirus species to identify distinct viral enzyme characteristics.

The cluster of differentiation 36 (CD36) is a membrane protein with broad physiological roles in health and disease, and its function is regulated in part by phosphorylation. Experimental evidence shows that phosphorylation of Thr92 reduces CD36 affinity for thrombospondin-1 (TSP-1), binding of which initiates antiangiogenic signaling, whereas phosphorylation of Ser237 decreases CD36-mediated fatty acid uptake, with implications for energy metabolism. However, the only available crystal structure of CD36 lacks phosphorylation, and the molecular mechanisms by which phosphorylation regulates CD36 function remain largely unknown. This study provides an atomically detailed computational characterization of CD36 in unphosphorylated and dual phosphorylated states, using molecular dynamics simulations with a total sampling time of 30 microseconds in combination with Markov state models. We present, to our knowledge, the first evidence of a cryptic pocket on CD36 surface that is formed by phosphorylation. This cryptic surface pocket and a loop spanning residues 121-131 form a high affinity binding site for TSP-1 derived ligands, shifting their binding away from the canonical site. We propose that this altered binding provides a molecular basis for the disruption of antiangiogenic signaling upon CD36 phosphorylation. Additionally, our data indicate that, phosphorylation increases helicity and compaction within the helix-loop region spanning residues 296-331, narrowing one of the entrances to the internal cavity and reducing its overall volume. These conformational changes provide a potential mechanistic explanation for the decrease in fatty acid uptake upon CD36 phosphorylation. Our findings provide structural insights that may inform the future design of CD36 modulators and emphasize the importance of targeting phosphorylation induced CD36 conformations in angiogenic and metabolic diseases.

Initiator-methionine excision is among the earliest protein modifications, yet its relationship to assembly geometry is unknown. Burial of the mature first residue was measured across 7,246 deposited human biological assemblies (22,291 chain-level observations; 1,191 proteins). Among 1,143 analyzable proteins, termini in MetAP-permissive penultimate-residue sequence classes were less often interface-engaged than termini in MetAP-nonpermissive classes (37.4% versus 47.4%; adjusted odds ratio 0.65, p = 7.2e-4). Curated processing annotations did not show a corresponding burial difference, and correlated residue properties preclude attributing the sequence-class association specifically to iMet removal. The analysis identified 264 interface-engaged MetAP-permissive candidates concentrated in cellular machines. In a fully recomputed conformer scan of deeply buried proteasome positions, modeled methionine accommodation was less favorable than at observed-methionine controls (median overlap -0.30 versus -1.12 angstrom, p = 0.0049), although most scoreable sites permitted a nonoverlapping placement. The census therefore reveals a graded structural constraint - not universal steric failure - and prioritizes complexes in which altered packing, assembly kinetics, lipidation or N-terminal methylation can be tested.

The folding thermodynamics of proteins are dominated by two opposing forces, the loss in backbone entropy and the packing of hydrophobic groups. The same forces are major contributors to the extension thermodynamics of elastic proteins with the distinction that both processes act in concert, favoring the higher chain and solvent entropy of a relaxed conformation. The relative entropic contributions specify the recoil mechanism; human elastin recoil is primarily driven by hydrophobic forces, whereas fly resilin has a rubber-like mechanism driven by backbone entropy. Despite the importance of elastic proteins to tissue biomechanics, few have been identified, let alone characterized to the same extent as elastin and resilin. We develop a thermodynamic framework that maps proteins by sequence-derived estimates of extension-induced backbone and solvent entropy changes. Putative elastic proteins are proposed and classified by recoil mechanism based on estimated thermodynamic features. Proteins that map to elastic regions are overrepresented by the skin proteome. The set of predicted elastic domains is further extended by incorporating sequence context embedded in protein language models. Protein domains with distinct thermodynamic recoil mechanisms cluster on the latent space manifold. Some of these domains are anticipated to have roles within molecular machines, expanding the scope of elastic protein function beyond mechanical materials like elastin and resilin.

We present the solution structures of HIV-1 proteins NC(p7)1-55 corresponding to the full-length NC(p7) and mature p6. The studies were carried in water and, to mimic the membrane, in micellar DPC (Dodecylphosphocholine) conditions. Our results unravel for the first time the structure adopted by the N-terminal amino acids of the free NC(p7)1-55, with the formation of a small helix spanning residues F6 to R10. Our NMR and Fluorescence Anisotropy data disclose an interaction between NC(p7)1-55 and p6 both in water and DPC, with respective Kd of 2.5mM and 370 mM at 23{degrees}C. The interaction is thus strengthened in lipidic conditions. Protein p6 stabilizes the N-terminus of NC(p7)1-55 while increasing at the same time the dynamic of the first zinc finger. Although the entire p6 sequence is involved in the interaction, we show that its C-terminal region is particularly sensitive to the presence of NC(p7)1-55, with a propensity of forming a a helix ranging from amino acids S111 to F116. This study brings experimental evidence of a direct protein-protein interaction between p6 and the N-terminal region of NC(p7)1-55. We further show that such interaction is readily accommodated within the NC(p15) framework and hypothesize that it may facilitate the selective assembly of assembly of the viral genomic RNA (gRNA) in the cell.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease characterised by progressive motor neuron degeneration. Mutations in the SOD1 gene represent the second most common genetic cause of ALS (ALS), and distinct SOD1 missense variants present with markedly different clinical profiles. A4V leads to an aggressive form of the disease (median survival [~]1y), H46R confers a mild, slowly progressive course and I113T exhibits an intermediate phenotype. The molecular basis by which these mutations produce divergent clinical outcomes remains poorly understood. We performed extensive classical molecular dynamics simulations of wild-type SOD1 and the three ALS-associated variants in the apo monomeric state to attempt to investigate the mechanisms behind such phenotypic differences. Structural stability, global compactness, and conformational flexibility, as well as analysis of collective motions between residues and estimation of free energy, were assessed. The H46R, A4V, and I113T variants exhibited distinct dynamic behaviours, highlighting differences in structural stability, local flexibility, and intramolecular interactions. These findings suggest that specific structural regions may contribute differently to protein dysfunction and could represent key elements for understanding the relationship between molecular dynamic properties and the differing clinical severity associated with these variants. Most strikingly, H46R exhibited exceptional structural stability across every analytical level, the lowest global deviation, most attenuated local flexibility, strongest internal dynamic coordination, and the deepest, most confined free energy basins of any system examined. This convergent multi-layered evidence of structural restraint provides a compelling mechanistic basis for the mild and slowly progressive clinical course of H46R ALS, suggesting that enhanced conformational rigidity, rather than bulk destabilisation, is the defining biophysical feature of this variant, and that its pathogenic mechanism operates through a route fundamentally decoupled from the aggregation-driven toxicity that characterises the more aggressive SOD1-ALS mutations.

Abstract Imaging-based spatially resolved transcriptomics (imSRT) technologies provide high-throughput molecular-resolution spatial characterization of genes within cells. Conventional analysis methods to identify cell-types and states in imSRT data rely on gene count matrices derived from tallying the number of mRNA molecules detected for each gene per segmented cell, thereby overlooking subcellular heterogeneity that can be useful in defining cell states. To take advantage of the molecular-resolution information in imSRT data and potentially identify cell-states based on subcellular heterogeneity, we developed STARIT (Spatial Transcriptomics As Rasterized Image Tensors). STARIT converts transcripts within segmented cells in imSRT data into an image-based tensor representation that can be combined with deep learning computer vision models for downstream analysis. Using simulated and real imSRT data, we demonstrate that STARIT distinguishes transcriptionally distinct cell-types and further separates cell states based on subcellular transcript localization, which conventional gene count analysis fails to capture. By providing a standardized framework to encode subcellular molecular information in imSRT data, STARIT will enable deeper insights into subcellular heterogeneity and enhance the identification and characterization of cell-types and states that are overlooked by gene count representations.

Meso2EM is a correlative light and electron microscopy workflow that transfers neurons selected from mesoscale functional images to targeted electron microscopy. We recorded Ca{superscript 2} signals from layer 2/3 neurons across a contiguous 3 x 3 mm cortical field in awake mice and reidentified a selected neuron after fixation and tangential sectioning. Lectin-labeled vascular architecture served as a shared landmark across in vivo two-photon imaging, confocal microscopy, laboratory micro-CT of resin-embedded tissue, and block-surface scanning electron microscopy, guiding focused-ion-beam scanning electron microscopy to the target cell body. The same progressive-targeting principle also supported serial ATUM-SEM reconstruction of an in vivo-tracked dendrite and serial transmission electron microscopy of optically selected dendrites from a patch-clamp-recorded Martinotti cell. Meso2EM therefore provides a practical route for preserving target identity across large changes in scale and specimen state while restricting electron-microscopy acquisition to a selected region.

Spatial-omics workflows enable molecular analysis within tissue spatial context. Despite the prognostic value of tissue stiffness, these approaches have not incorporated direct, mechanical measurements. This omission reflects several challenges, including sample requirements, low throughput, specialized equipment, and complex data registration. Here, we introduce mechanoMaST (mechanics mapped to spatial transcriptomics), the first workflow to combine absolute mechanical measurements with spatial-omics. It pairs atomic force microscopy-based nanoindentation stiffness maps with spatial transcriptomics maps from adjacent tissue cryosections. The two modalities are then computationally co-registered to enable direct spatial correlation at 100 um resolution, with mapping accuracy quantified through error propagation, providing ground-truth mechanical data directly linked to spatial gene expression. We demonstrate mechanoMaST in human colorectal cancer liver metastasis, generating a spatial resource from 10 patients and revealing a four-gene stiffness signature. mechanoMaST is readily adaptable to other tissues across development and disease, and extendable to additional spatial-omics modalities in adjacent sections.

Deep learning-based protein structure prediction methods that leverage evolutionary information from multiple sequence alignments (MSAs), exemplified by AlphaFold2, have achieved remarkable accuracy. However, existing methods still struggle to predict challenging proteins, particularly those with novel folds or limited evolutionary information, and to recover alternative conformational states. Here we show that structure prediction models trained under different MSA-depth distributions corresponding to different levels of evolutionary information exhibit complementary generalization behaviors, and that a model trained on a mixture of these distributions can combine their complementary generalization strengths. Building on this insight, we developed ProtMonomer, a deep learning framework trained on MSA-depth distributions representing a broad range of evolutionary information levels to improve structure prediction. Across benchmarks comprising CASP15 targets, non-redundant experimentally determined structures, orphan proteins, and short peptides, ProtMonomer performed comparably to or better than leading methods, including AlphaFold2 and AlphaFold3, with particularly strong performance on challenging targets. For fold-switching proteins, ProtMonomer also recovered alternative conformational states more accurately than AlphaFold2 and AlphaFold3 across diverse homologous sequence sampling strategies. In addition to improving predictive accuracy, ProtMonomer substantially reduced inference cost through an efficient architecture, enabling high-throughput applications. Together, these findings provide insights into the generalization of evolution-informed structure prediction models and support ProtMonomer as an accurate and efficient framework for protein structure prediction.

Surface phenotyping underpins plant science, preclinical animal research and entomology, yet across all three the measurement is almost always a photograph, which records a projection and not the surface itself. Here we present the Gentschinator3000, an open structured-light platform that brings high-end metric surface measurement within reach of laboratories with no optics expertise, combining documented open hardware, open reconstruction software and analysis workflows for under 4000 Euro in components. It resolves a planar reference to 45 m local flatness, registers full rotations to a loop closure of 156 m, and performs stably across acquisition ranges that we define. Applying one workflow to a leaf before and after desiccation, to murine anatomy and to a spread lepidopteran, we find that projection underestimates surface area by 11 to 41 %. That error grows with the condition under study, with the evaluation scale and with the direction of view, so it can confound phenotype comparisons dramatically. In murine limbs a 15-degree change of viewing direction shifts a projected inter-segment angle by up to 23.2 degrees, while the three-dimensional angle does not move. Projection geometry can therefore contribute as much to a measured phenotype as the biology it is meant to quantify.

Linker chemistry and conformation are central determinants of PROTAC activity, shaping ternary-complex geometry, cooperativity, target-lysine presentation and cellular permeability. Existing linker generators often lack explicit control over linker flexibility, require predefined attachment sites and linker lengths, or produce structures that demand substantial geometric correction, limiting their utility in practical PROTAC design. Here we introduce FlexiTAC, a Bayesian flow network that jointly generates linker atom types and coordinates from the warhead and E3-ligase-ligand contexts. We also assemble PROTAC-3D, a quality-controlled collection of 63,554 component-resolved PROTAC structures for model training, and PROTAC-Bench, which covers molecular quality, fragment preservation, geometric fidelity, conformational stability, fragment awareness, rediscovery and sampling efficiency. Compared to the best 3D baseline models, FlexiTAC improves validity by 12.0-12.7% and achieves the highest PoseBusters pass rate of 79.5%-80.0%. A differentiable guidance module shifted generated linkers along a conformational ensemble-derived rigidity axis without retraining the generator. In silico case studies further show that the model can accept crystal-derived, redocked or predicted structural inputs. Together, FlexiTAC, PROTAC-3D and PROTAC-Bench establish an integrated and reproducible framework for data-driven PROTAC linker design, combining controllable structure-conditioned generation with standardized training data and evaluation protocols. This framework expands the linker chemical and conformational space accessible to computational exploration, provides a foundation for future method development and enables the systematic generation of structure-conditioned linker designs with tunable conformational flexibility.

Label-free Raman imaging of living cells is photon limited: at exposures compatible with cellular dynamics, single-pixel spectra carry about one count per channel on a dominant smooth background. We present an unmixing framework in which the decoder of a physics-constrained autoencoder is restricted to a data-driven spectroscopic dictionary: band centers,widths, and pseudo-Voigt shapes are measured from the dataset and fixed, and the network learns only nonnegative band amplitudes, a smooth B-spline background, and a per-pixel gain.First, on slit-scanning images of HeLa cells (532 nm) the dictionary yields spike-free component spectra that read as band tables, including a resonance-enhanced cytochrome-c-associated component matching literature spectra, and the most stable decomposition against the component number. Second, the dictionary and initialization calibrated at 1 s exposure perline transfer to 100 ms per line (12 s sweeps): cytochrome-c spectral identity survives a single sweep (correlation 0.92) while its map remains photon limited; the dictionary provides spectral physicality, and the transferred initialization prevents a structural collapse that global map correlations miss; in a measurement-derived phantom the dictionary estimator holds thecytochrome-c spectrum to 17-19{degrees} spectral angle at 100 ms, where classical factorizations and free decoders lose it (55-64{degrees}). Estimation on the count-equivalent detector output uses a calibrated shifted-Poisson quasi-likelihood. Third, evaluation must be time matched:correlation against a separately acquired reference saturates through slow specimen drift and acquisition mismatch rather than photon noise, and the self-consistency of learned denoisers is inflated by shared bias; time-matched self-consistency and independent cross-checks areproposed.

Transcription factor Ets2 coalesce with the NF-{kappa}B pathway to regulate gene expression in specific signaling contexts. IKK2/{beta}-mediated phosphorylation events critically regulate the NF-{kappa}B pathway. However, any link between Ets2 and IKK2 remains elusive. Here we report Ets2 as a direct substrate of IKK2. In-vitro kinase assays using deletion constructs, high resolution MS-MS and site directed mutagenesis identified S295 as a prominent phosphorylation site distal to the DNA binding domain, substitution of which to phosphor-mimetic Glutamate triggers further phosphorylation of Ets2. MD simulations clearly indicate conformational constriction of the otherwise disordered N-terminal region and inhibition of DNA binding activity upon phosphorylation, which was further confirmed by Electrophoretic mobility shift assays. Our results uncover a phosphoregulatory connection between Ets2 and IKK2.

Recent advances in protein structure prediction, exemplified by AlphaFold, have largely addressed the determination of static structures, one aspect of the protein folding problem. However, predicting folding pathways, by which proteins reach their native states, remains a significant challenge. Here, we present PathFold, a deep learning framework that predicts protein folding pathways directly from sequence information. PathFold leverages an AlphaFold-based module to extract structural information from the sequence and generates a progressive folding trajectory from an extended conformation using a diffusion model. By modeling the full trajectory, it enables prediction of folding intermediates and transition pathways, analogous to those observed in steered molecular dynamics (SMD) simulations. The predicted pathways reveal well-defined intermediates and sequential folding events, and show agreement with experimental folding data, including measured {Phi}-values.

Despite the well-known health risks of neutron exposures, key gaps remain in understanding neutron-induced molecular responses and identifying reliable biodosimetric markers that distinguish neutrons from photon exposure. We provide the first genome-wide analysis of the human blood transcriptional response to an accelerator-derived fission-like spectrum of neutrons versus photons, evaluating transcriptomic relative biological effectiveness (RBE) and radiation quality-discriminating gene signatures. Whole blood from healthy donors was irradiated ex vivo with X-rays (140 kV, 0-4 Gy, n = 3) or neutrons (0.1-8 MeV, 0-1 Gy, n = 2), incubated for 6 h or 24 h, and processed for RNA sequencing from peripheral blood mononuclear cells (PBMCs). Neutrons were markedly more potent than X-rays at inducing differentially expressed genes (DEGs) at equal doses, showing a peak response 6 h post-irradiation followed by a decline. In contrast, X-rays caused a continuous increase in DEGs up to 24 h (neutrons vs. X-rays at 1 Gy: 1,449 vs. 121 DEGs at 6 h; 996 vs. 621 DEGs at 24 h). A universal p53-centered 34-gene signature, including FDXR, EDA2R, GADD45A, and ZMAT3, showed highly monotonic dose responses (Spearman correlation coefficient {approx} 1) across donors, radiation qualities, and timepoints. Additionally, difference-in-differences analysis identified radiation quality-discriminating genes only at 6 h, with transcriptional convergence observed by 24 h, suggesting a very narrow time window for biodosimetric differentiation. We identified a neutron-specific gene signature driven by cGAS-STING-NF-{kappa}B signaling (RELB, NFKB1, C3, MALAT1) and suppression of B-cell and myeloid identity genes (IGHD, TCL1A, CLEC7A, TLR2), defining a biologically coherent neutron quality index with distinct immunomodulatory effects. For the first time, we assessed neutron RBEs at the gene, pathway, and global transcriptomic levels in a human blood model, reporting a global transcriptomic neutron RBE of 1.30 (95% CI: 1.14-1.49) at 6 h and 1.21 (95% CI: 1.14-1.28) at 24 h, providing a valuable basis for biodosimetry in mixed-field exposure scenarios. Our findings advance the mechanistic understanding of neutron radiation responses and support the development of biodosimetric approaches for mixed-field exposure scenarios.